Transrelativistic pair plasmas in AGN jets 
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ABSTRACT 

Models of relativistic jets filled with ultrarelativistic pair plasma are very successful 
in explaining the broadband radiation of 7-ray blazars. Assuming that the initial 
injection and cooling of ultrarelativistic pair plasma in an AGN jet has occurred, 
producing the observed high-energy 7-ray radiation, we investigate the further 
evolution of the pair plasma as it continues to move out from the central engine. The 
effects of thermalization and reacceleration, the emission of pair bremsstrahlung and 
annihilation radiation and the bulk Compton process, and the possible application to 
MeV blazars are discussed. A model calculation to the special case of PKS 0208-512 is 
presented. 

Subject headings: acceleration of particles — plasmas — radiation mechanisms: 
non-thermal — radiation mechanisms: thermal — radiative transfer — galaxies: jets 
— gamma-rays: theory 

1. Introduction 

The detection of high-energy 7-ray emission from more than 60 blazars with EGRET is a 
challenge and at the same time a constraint of fundamental importance for emission models (von 
Montigny et al. 1995). A large fraction of these blazars exhibits variability at 7-ray energies on 
time scales of days to months (Mukherjee et al. 1997). The optical counterparts of the majority of 
EGRET detected AGN are known as BL Lacertae objects and optically violent variable quasars 
(OVV). At radio wavelengths, all blazars can be recognized as bright, compact sources with a 
flat synchrotron spectrum emanating from outflowing plasma jets that are nearly aligned with 
our line-of-sight. Relativistic beaming is required in the objects in view of the luminosity and 
variability time scales (Dermer and Gehrels 1995), in accord with VLBI observations indicating 
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that superluminal motion is a common feature in this class of AGN (e.g. Wehrle et al. 1994, Pohl 
et al. 1995, Barthel et al. 1995, Krichbaum et al. 1995). 

The strongest EGRET blazar detections can be characterized by a single power-law spectrum 
with differential photon spectral indices between a = 1.5 and a = 2.7 (Thompson et al. 1995). 
For individual sources the spectral index is correlated with the flux level, and there may also 
be deviations from the power-law behaviour both below 70 MeV and above a few GeV (Pohl et 
al. 1997). The combined OSSE/COMPTEL/EGRET measurements generally indicate spectral 
breaks at a few MeV (Williams et al. 1995, McNaron-Brown et al. 1995) in the sense that the 
spectra below 1 MeV are harder than in the EGRET range. 

At medium 7-ray energies observable by COMPTEL, PKS 0208-512 has been identified as 
an AGN with flaring properties at MeV energies (Blom et al. 1995). There is now evidence that 
it may belong to a class of 'MeV-blazars' that are occasionally exceptionally bright MeV sources 
(see Bloemen et al. 1995). Since in these objects the bright emission appears to be confined to a 
relatively narrow energy range, the discussion has focused on models involving a broad blue-shifted 
e + /e~ annihilation line that is Doppler boosted in a relativistic jet (Roland and Hermsen 1995). 
If this interpretation applies, MeV-blazars provide a unique tool to study astrophysical particle 
beams (e.g. Schlickeiser 1996). Furthermore, MeV-blazars may confuse the analysis of galactic 
sources of 7-ray line emission (Pohl 1996). 

In an earlier paper (Bottcher, Mause & Schlickeiser 1997; hereafter BMS) we have investigated 
the temporal evolution of ultrarelativistic pair plasmas in jets of quasars and BL-Lac objects and 
have demonstrated that their broadband spectra can well be explained as the resulting synchrotron 
and inverse-Compton radiation from a cooling ultrarelativistic nonthermal pair plasma in a 
relativistic jet. A decisive parameter for the evolution of single ultrarelativistic plasma components 
inside an AGN jet is the density of pairs injected into the jet. Broadband fits, covering the radio 
to 7-ray regime of the electromagnetic spectrum, to blazars generally require particle densities 
of order n e < 10 3 cm -3 . Nevertheless, fits to different objects suggest that the value of the pair 
density in relativistic jets ejected by active galactic nuclei varies over several orders of magnitude. 
Jets of very high density (n > 10 5 cm -3 ) can also produce the broadband spectra at least of 7-ray 
active flat-spectrum radio quasars whose bolometric luminosity is clearly dominated by the 7-ray 
emission. 

In this paper, we argue that after an initial phase of rapid cooling, governed by synchrotron 
and inverse-Compton energy losses, the relativistic pair plasma inside such components is likely to 
attain a quasi-thermal distribution. If the jet remains well collimated and cooling (e. g. through 
adiabatic losses) remains very efficient through the transrelativistic phase, pair bremsstrahlung 
and pair annihilation become efficient. The resulting radiation spectrum peaks around several 
MeV, which has been suggested previously to be responsible for the observed MeV bump in 
MeV blazars (Henri et al. 1993, Roland & Hermsen 1995, Bottcher & Schlickeiser 1996). 
However, while generally the high-energy radiation from a few ultrarelativistic jet components, 
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ejected over the typical EGRET variability time scale, is appropriate to model EGRET spectra, 
the pair annihilation and bremsstrahlung radiation from only a few components after cooling 
and thermalization is usually too weak to explain the MeV blazar phenomenon. Therefore, a 
quasi-continuous supply of mildly relativistic pair plasma into the jet is required in order to fit the 
MeV bump in MeV blazars with pair annihilation radiation. 

Alternatively, if very close to the central accretion disk the inverse-Compton cooling rate is 
balanced by reacceleration by hydromagnetic turbulences the pair plasma in the jet will thermalize 
at relativistic temperatures and will continuously Compton upscatter external radiation from the 
accretion disk very efficiently, producing another bump at keV — MeV energies by the so-called 
bulk Compton process (Sikora et al. 1997). 

We know from the EGRET data of many blazars that the GeV emission is not that of a single 
injection of particles, but is more likely to result from more or less regularly repeating injection 
events with varying energy input into the ejected particles. A second important point is the 
correlation between flux and spectral index in the EGRET range. The spectral hardening during 
outbursts indicates a more efficient acceleration of particles. Even for highly variable sources like 
PKS 0208-512, PKS 0528+134 (Collmar et al. 1997) and others the low-energy 7-ray continuum 
varies with small amplitudes and not in phase with the variability in the EGRET range. To be 
more realistic we may thus either assume that injections occur in regular time intervals, so that 
they would not influence each other in their evolution, or that injection occurs in a quasi-steady 
manner, such that relativistic pair plasma is continuously injected. 

For ease of both computing and exposition we will not consider the effect of the finite light 
travel time within the volume occupied by pair plasma, which limits our predictions of variability 
to time scales longer than the light travel time through individual plasma blobs. In section 2, we 
give a short overview of the different elementary processes which play a role for the evolution of a 
relativistic pair plasma component of an AGN jet after the initial phase of rapid cooling and which 
are usually not considered in models of ultrarelativistic jets. In section 3, we describe numerical 
simulations and an analytical approximation, applicable under special conditions, to follow the 
evolution of the pair plasma through the transrelativistic phase. We find that, depending basically 
on the particle density and the magnetic field, there are two different ways how a quasi-thermal 
distribution can be established. The relevant physics of these quasi-thermal plasmas is described 
in section 4. Section 5 contains a model calculation for the typical MeV blazar PKS 0208-512. We 
summarize in section 6. 

2. The elementary processes 

The momentum distribution of pairs inside a relativistic jet component is assumed to be 
isotropic and (for sake of simplicity) homogeneous. Under circumstances which will be specified 
below, the time evolution of such a pair plasma can be approached by a Fokker-Planck equation 
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regarding only time and particle energy as variables: 
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Here, ^ is the single-particle energy loss due to all the relevant elementary processes and is 
the energy dispersion rate. The term on the right-hand sight of Eq. (1) is the source function. We 
restrict ourselves to the case where the plasma component is optically thin to 7-7 pair production 
(this process might well have led to re-injection of pairs during the ultrarelativistic phase but 
becomes irrelevant in the transrelativistic regime), and consider the catastrophic pair annihilation 
losses and the dilution of the blob due to expansion. 

A detailed description of the calculation of the energy loss/gain rates due to inverse-Compton 
scattering of accretion disk radiation, synchrotron emission and the SSC process can be found in 
BMS. At the final states of the simulations carried out there, the particles have cooled down to 
energies such that all Compton scattering events may be treated in the Thomson regime. This 
guarantees that the use of the Fokker-Planck equation to describe the effect of inverse-Compton 
scattering is a good approach, which requires that the change in particle energy is A7 <C 7 — 1 for 
each scattering event. Since A7,< e (j 2 — 1) = e( 7 + 1) (7 — 1) and 67CI, this condition is fulfilled 
in the case of Thomson scattering. The external photon source is assumed to be dominated by 
the accretion disk. The full angle-dependence of the accretion-disk photon field is included in our 
calculation. 

Modelling of blazar spectra indicates that the magnetic fields in many AGN jets can be 
significantly lower than the equipartition value. Therefore, the gas pressure generally exceeds the 
magnetic pressure, leading to a conical jet and freely expanding blobs. The inferred low magnetic 
fields (typically B < 0.1 - 1 G) in some objects, together with the fact that the magnetic field 
generally declines outward, imply that synchrotron radiation (and therefore also reacceleration by 
synchrotron-self absorption) is of minor importance in the evolutionary phase with which we are 
dealing in this paper. 

In ultrarelativistic jets, the effects of elastic (M0ller and Bhabha) and inelastic scattering 
(pair bremsstrahlung emission), adiabatic losses, pair annihilation and stochastic reacceleration 
may be neglected compared to the dominant influence of synchrotron and inverse-Compton losses 
(BMS). Being interested in the details of the process of cooling down to mildly-relativistic energies, 
we now have to consider them and therefore discuss their influence in the following subsections. 
In this phase of the jet evolution, also the expansion of the jet will have a significant influence via 
adiabatic cooling and the dilution of the particle densities, affecting elastic and inelastic scattering 
and pair annihilation. 
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2.1. Elastic scattering 

The energy loss/gain rate due to M0ller scattering of electrons (positrons) of energy 7$ off 
an isotropic distribution n±( r )±) = 47r.ZV±/?±7± /±(7±) of positrons (electrons) in the general case 
has been calculated by Dernier (1985). Here, we may restrict our considerations to the relativistic 
case in which the elementary single-particle energy loss is dominated by the term involving the 
Coulomb logarithm. 

As was suggested by Nayakshin & Melia (1998), in this case the single-particle energy 
gain/loss rate reduces to a one-dimensional integral: 
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A = lnA + -ln(2e) « lnA + 0.8466. (4) 

The Coulomb logarithm In A is only weakly dependent on the electron and positron energies. In 
our numerical calculations, we use the constant value A = 20. The difference between M0ller 
(electron-electron) and Bhabha (electron-positron) scattering only occurs in the terms neglected 
in the derivation of eq. (3) and is negligible in the case of relativistic pairs. The energy-exchange 
rate scales linearly with particle density and, for energies smaller than the average particle energy 
of the distribution, 7 < (7), it declines roughly as dj/dt oc 7 _1 . It changes the sign at 7 = (7) 
and approaches a constant, negative value for 7 S> (7). 

In order to evaluate the energy dispersion due to elastic (M0ller and Bhabha) scattering, 
we use the technique outlined by Dermer (1985), together with the idea of Nayakshin & Melia 
(1998) to consider scattering off monoenergetic distributions (for which the dispersion rate can 
be calculated analytically) and then averaging over the distributions of the background particles. 
In contrast to the technique of Dermer (1985), this requires only one numerical integration. The 
resulting expressions are given in Appendix A.l. The energy-dispersion rate also scales linearly 
with particle density and depends only weakly on the particle energy in the range I7 — (7)) <C (7). 
The Fokker-Planck equation (1) is only valid if the energy distribution is narrow in the above 
sense, because otherwise elastic scattering events between particles of very high and very low 
energies would lead to a considerable energy gain A7 >, 7 of the low-energetic particle. 
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2.2. Bremsstrahlung emission 



The general treatment of the pair bremsstrahlung process is very cumbersome and has been 
investigated in several papers by Haug, posing special interest on pair bremsstrahlung emission in 
thermal pair plasmas (Haug 1975, 1985a, 1985b). The importance of this radiation mechanism to 
the hard X-ray and 7-ray emission of relativistic jets in AGNs has been demonstrated by Bottcher 
& Schlickciser (1995). 

For the pair distributions considered in this paper, we may restrict ourselves to the 
ultrarelativistic limit using the differential pair bremsstrahlung cross section derived by Baler, 
Fadin &; Khoze (1967). The bremsstrahlung energy-loss rate is calculated in Appendix A. 2. 
Inspecting the function K given there, we find that the final-state averaged energy-loss of a test 
particle emitting bremsstrahlung photons is typically of the order m e c 2 and thus, for relativistic 
particles the necessary condition A7 <C 7 — 1 is fulfilled. The energy-loss rate has the same 
linear dependence on particle density as for elastic scattering and depends on particle energy as 
—dj/dt oc 7 L1 . It is only weakly dependent on the shape of the relativistic pair distribution. 



For the parameter range under consideration in this paper, we expect the jet to expand freely. 
This yields a declining particle density as 



The expansion also leads to adiabatic losses due to the invariance of the magnetic flux through 
the particle orbits, Br\ = const., where ri is the Larmor radius of a particle. This provides an 
additional energy-loss term, 7 = B/(2B), depending on the z dependence of the magnetic 
field. In the case of a perfectly isotropized magnetic field, its evolution is given by B(z) oc z~ 2 , 
while a purely transversal magnetic field evolves as B(z) oc z^ 1 (Blandford &; Rees 1974, Blandford 
& Konigl 1979). Since the degree of isotropization of the magnetic field is not well known, we 
parametrize the magnetic field dependence as 



2.3. Effects of jet expansion 
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with 1 < b < 2. This yields 



2.4. Stochastic acceleration by Alfven waves 



In order to estimate the reacceleration and energy dispersion rate due to wave-particle 
interaction we will restrict our calculation to parallel, transverse plasma waves, i. e. k || Bo and 
E_LBq, in a pure pair plasma. Here, k is the wave vector, E is the polarization vector of the 
plasma wave and Bo is the background magnetic field. 

According to the scenario of the formation of jet components, described in detail in BMS, 
the role of protons inside the jet may be neglected. In this scenario, the pairs are created as 
secondary particles in photo-pair and photo-pion production by relativistic protons accelerated in 
the accretion disk magnetosphere and will greatly outnumber the primary particles and carry the 
dominant fraction of the total energy transferred to particles, if the duty cycle of 7-ray emission 
is much less than unity. This pair population becomes unstable with respect to the excitation of 
various electromagnetic and electrostatic waves, resulting in an explosive event which ejects the 
ultrarelativistic pair plasma along an existing jet structure. This implies that the mechanism for 
acceleration of the protons during the quiescent phase by stochastic processes, ultimately leading 
to the production of the unstable pair plasma population, is of completely different nature than 
the mechanism which injects the pair plasma into the jet. The magnetic field and the level of 
turbulence during the latter phase, in turn, are in general drastically different from their values 
along the jet structure. 

It is well-known that relativistic particles can only interact efficiently with the Alfven part 
of the plasma wave spectrum. For this reason, we neglect interaction with electron (positron) 
cyclotron waves. The relevant Fokker-Planck coefficients for the interaction of electrons with 
Alfven wave turbulences propagating through an electron-proton plasma have been investigated 
in detail by Schlickeiser (1989). Here, we generalize his calculation to the case of Alfven waves in 
a pair plasma. We assume that the Alfven waves have wave numbers k between 

k ■ -— and k -— (9) 

">mm — n a±±u "'max — \J ) 

Rb v a 
where Rb is the radius of the jet component, Qq is the nonrelativistic gyrofrequency of 
electrons/positrons, and v a = Bq v 7 (7) / V 87rm e n e is the Alfven velocity in a relativistic pair 
plasma. 

Eq. (9) certainly underestimates the minimum wave number since we assume the magnetic 
field to be ordered on smaller scales than Rb- It is especially the highest-energetic particles 
which interact with these longest waves. But for ultrarelativistic particles, synchrotron and 
inverse-Compton cooling is much more important than wave-particle interactions. Thus, the 
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inaccuracy of Eq. (9) does not influence our final results. We assume that the spectral energy 
content in plasma waves is distributed according to a power-law, 



with 



I{k)=hk~ <l (10) 
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Here, 5B is the amplitude of the magnetic field fluctuation due to the plasma waves which is 
typically of order <, lO^ 1 Bq. Now, following the calculations of Schlickeiser (1989) (accounting 
for the differences of the dispersion relation between an e~-p- and an e~-e + -plasma), we find the 
acceleration rate due to Alfven waves, 

and the energy dispersion rate 
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2.5. Energy dispersion due to inverse- Compton scattering 

Being a stochastic process, also inverse-Compton scattering leads to energy dispersion. The 
respective energy-dispersion rates for Compton scattering in the Thomson regime are calculated 
in Appendix A. 3. 

For relativistic particles, the energy-dispersion rate due to inverse-Compton scattering is 
d(A r y) 2 /dt oc 7 4 and linearly dependent on J °° de e 2 n p h(e), i. e. the photon densities and the 
average of the square of the photon energies. Their extremely strong dependence on particle 
energy (oc j 2 and oc 7 4 , respectively) implies that the energy-exchange and dispersion rates due to 
inverse-Compton scattering will usually dominate the total rates for very high particle energies. 

Which process is dominant at low particle energies, depends critically on the particle density 
and the strength of magnetohydrodynamic turbulences. In the case of high densities, elastic 
scattering will govern the evolution of these particles, while in the case of lower density and strong 
plasma wave turbulences, wave-particle interactions (reacceleration) and adiabatic losses will 
dominate their behavior. 

Figs. 1 and 2 show the energy-exchange and energy-dispersion rates for the different processes 
for a plasma of density 3 • 10 6 cm~ 3 , having a narrow distribution around (7) « 30. A magnetic 
field of B = 0.1 G is assumed and the blob is located at z = 10 -3 pc above an accretion disk 
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with total luminosity L = 10 45 ergs -1 and moves outward with bulk Lorentz factor V = 15. The 
magnetic-field is assumed to be isotropized (6 = 2). These parameter values are similar to those 
which result from the cooling of an ultrarelativistic pair plasma producing an X-ray and 7-ray 
spectrum consistent with the observed spectrum of PKS 0208-512. 



2.6. Pair annihilation losses 

A detailed discussion of the effects of pair annihilation on a relativistic pair plasma in 
AGN jets can be found in Bottcher & Schlickeiser (1996). In order to calculate the total pair 
annihilation rate, we use the expression given by Svensson (1982) which, for completeness, is 
quoted in Appendix A. 4. For relativistic particles, it declines roughly as n±(7±) oc (7±(7:p)) _1 , 
where (7^) is the average Lorentz factor of electrons or positrons, respectively, and becomes only 
important in the case of very high particle densities. 



3. Simulations of the jet evolution 

Since several energy loss and dispersion rates depend on the present particle distributions, Eq. 
(1) in its general form represents a set of coupled integro-differential equations which, in general, 
can not be be solved analytically. For the purpose of numerical integration, we used an explicit 
code of finite differencing. Because the numerical simulations are extremely time-consuming 
we tried to find simplified approximative expressions for the coefficients in Eq. (1) in a way 
that the problem could be solved analytically. This is possible if both the cooling rate and the 
energy-dispersion rate are dominated by inverse-Compton scattering for all occupied energy states, 
which is the case if the pair distributions still have a relativistic low-energy cutoff at 7 ^ 10. Then 
pair annihilation is negligible and provided that the cooling time scale is much shorter than the 
time scale of dilution of the pair plasma due to expansion, the source term in Eq. (1) may be 
neglected. Furthermore, the energy-loss and dispersion coefficients are dominated by synchrotron 
emission and inverse-Compton scattering and can very well be approximated by 

f ~-A lX (tW, (14) 

where A{ are constants and both Fokker-Planck coefficients have the same explicit time 
dependence x(i). For the resulting differential equation we found an approximative analytical 
solution which is derived in Appendix B. Given a known distribution 71(70, 0) at t = 0, we can 
forward it in time to t = ti by 
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Since the explicit time dependence x(t) of the coefficients, being governed by the evolution of the 
dominant soft photon field, is difficult to separate in realistic situations, we choose time steps t\ 
in a way that we may assume the coefficients A\x{t) and Azx{t) to be constant within this time 
interval, i. e. we approximate x(t) by a st e P function. We use the solution (16) as long as the 
exact energy loss and dispersion rates of the lowest-energetic particles do not deviate more than 
10 % from the representations (14) and (15). Comparison with numerical integrations of Eq. 
(1) showed good agreement. After this phase we continue with numerical integrations of the full 
Fokker-Planck equation. 

It is obvious that we cannot extend the simulations using the code described above to 
subrelativistic pair temperatures since especially the bremsstrahlung cross section which we use 
is only valid in the ultrarelativistic limit. However, under a wide range of assumptions, the pair 
plasma attains a quasi-thermal distribution whose temperature depends basically on the energy 
density of the external radiation field, the pair density, the magnetic field and the amplitude of 
plasma wave turbulences. 

Figs. 3 — 6 show two examples of this evolution as results of the numerical scheme described 
above. These examples illustrate the general tendency towards the establishment of a quasi- 
thermal distribution inside the jet component. At high particle energies the real quasi-equilibrium 
distributions are truncated, if radiative energy-loss processes dominate over adiabatic losses (Fig. 
4). If adiabatic losses dominate, the high-energy end of the particle spectrum tends to be slightly 
harder than a thermal distribution (Fig. 6). In both examples, the initial conditions were chosen 
in a way that the 7-ray compactness is low, which results in a time-averaged power-law 7-ray 
spectrum of photon number index a 7 = (s + 2)/2 due to inverse-Compton scattering of soft 
photons, where s is the spectral index of the electron energy distribution. 

In the case of extremely high particle densities (n > 10 8 cm -3 ), quasi-thermalization is a 
consequence of elastic scattering becoming the dominant energy exchange process. One example 
for this fact is illustrated in Figs. 3 and 4. Fig. 4 demonstrates that for such high densities, the 
truncation of the distribution functions with respect to a thermal distribution at high energies is 
of minor importance. 



For the case of lower densities, Schlickeiser (1985) demonstrated that if the spectral index q of 
the Alfven wave turbulence equals 2, the combined action of stochastic acceleration and radiative 



losses also establishes a quasi-thermal distribution. Its temperature is given by 



9 = ^/l + (3 + a)2p§ - 1 (18) 

where p c m e c is the particle momentum for which the net energy loss or gain, respectively, 
vanishes, and a = cui/a^ is the ratio of the coefficients governing the stochastic acceleration and 
energy-dispersion rates via dp/dt = u\p and d(Ap) 2 /dt = a2P 2 - The equilibrium temperature in 
Eq. (18) obviously depends very sensitively on the efficiency of radiative cooling and can become 
quite high if the external photon energy density is low, because the energy dependence of adiabatic 
(oc 7) and bremsstrahlung (oc 7 1 ' 1 ) energy losses is much weaker than that of inverse-Compton 
losses, leading to a very high equilibrium Lorentz factor (obviously, elastic-scattering const.] 
cannot balance acceleration at all). 

In our simulations, we find a similar behavior also for the case of Kolmogorov turbulence, 
i. e. for q = 5/3. Figs. 5 and 6 illustrate one example of this evolution. In this case, however, 
the deviation from a maxwellian distribution at high energies is more pronounced than in the 
case of a dense plasma. Assuming that inverse-Compton scattering of accretion disk radiation 
is the dominant cooling mechanism, which is balanced by stochastic reacceleration, and that 
quasi-thermalization occurs at relativistic temperatures, 9 > 1, we find for the quasi-equilibrium 
temperature 



e eq « 3 • 10 5 Bl nj 3 6^ -gL, (19) 

-^46 

where -Bo is the magnetic field in Gauss, -R15 is the blob radius in units of 10 15 cm, n e is in cm~ 3 , 
(SB/Bq) = 0.1(5_i, z pc = z/(lpc), and L46 = Le>/(10 46 ergs -1 ). In order to derive this estimate, 
we have used Eqs. (32), (37) and (38) in the limit 9 > 1. 

The X-ray and 7-ray emission of a quasi-thermal, mildly relativistic plasma is dominated by 
its pair annihilation and pair bremsstrahlung radiation. For pair temperatures < 3 the radiative 
output near the high-energy cut-off of the photon spectra is dominated by pair annihilation 
radiation which is insensitive to the pair distribution near its high-energy tail. The X-ray 
spectrum, to which only pair bremsstrahlung radiation contributes significantly and which has a 
photon number spectral index of a ps 1.1, is equally insensitive to the detailed shape of the pair 
distribution at high energies. For this reason, the observeable radiation from a mildly relativistic 
quasi-thermal pair plasma may well be approximated by the respective emission of a thermal 
plasma whose distribution coincides with the exact one up to the mean particle energy. 

This situation changes if the plasma maintains a highly relativistic temperature (0^5). 
In this case, the luminosities in synchrotron, SSC, and/or external Compton radiation may well 
dominate over pair bremsstrahlung and pair annihilation radiation, and the exact pair distribution 
has to be used in order to calculate the emitted X-ray and 7-ray spectrum. 
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The pair annihilation rate in a thermal plasma can be aproximated by 

dn 2 it c r 2 n 2 f 1 G 
~dt ~ (1 + 9) \l + 66 + 0.25 + ln(l + 1.15456) 
(Svensson 1982). The bremsstrahlung emission of a thermal pair plasma can be calculated as 

where g(k, 6) are the respective gaunt factors for e-e and e + -e~ bremsstrahlung for which we 
use the approximations for transrelativistic pair plasmas found by Skibo et al. (1995). The pair 
annihilation spectrum emitted by a thermal plasma is 

(Dermer 1984) where 
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(23) 



The bulk Compton process is calculated in the Thomson regime, including the full 
angle-dependence of the external radiation field (see BMS for details). 

We will now derive an estimate for the density required for elastic scattering to dominate the 
effect of stochastic reacceleration. To this aim, we compare the energy-loss due to elastic scattering 
for particles of energy 7 3> (7) to the respective reacceleration rate. Using Eq. (2), we find 



-6-lO^s- 1 (24) 

« t /Mo,7»( 7 > (7/ 

where 77,4 = n e /(10 4 cm -3 ). The stochastic acceleration rate according to Eq. (12) is 

|) a - 5 ■ 10- ^ R- 5 ^ (7)- 7 2/3 s~\ (25) 
This yields the condition 



n e > 3 • 10 7 B 7 /6 5-! R^ 1/3 7 1/3 cm -3 (26) 

for elastic scattering to dominate the energy loss at particle energy 7 » (7). If this condition is 
fulfilled for particles at the high-energy tail of the pair distribution, reacceleration is inefficient 
as compared to elastic scattering. In this case, radiative and adiabatic cooling will prevent the 
particle distributions from attaining a stationary state. 
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4. Evolution of a thermal plasma 

As we have shown in Section 3, in the case of very high particle densities (n > 10 8 cm -3 ), 
the pair plasma of a relativistic jet component approaches a quasi-thermal distribution due to the 
dominant action of elastic scattering. The evolution of a thermal distribution of particles, subject 
to pair annihilation, bremsstrahlung, inverse-Compton, synchrotron, SSC, and adiabatic losses, 
can be followed using the expressions for the respective cooling rates of thermal pair plasmas. 

For this purpose, we assume that the thermalization timescale (due to elastic scattering) is 
much shorter than the energy loss timescale. In this case, M0ller and Bhabha scattering will 
always establish a thermal distribution, 



dn{ 1 )= e--,h. (27) 
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where = kB \ is the dimensionless temperature and Ko is the modified Bessel function of second 
kind of order 2. In order to determine the respective cooling rates in a thermal plasma, we note 
first that 



K (- 

p 2 ) = (m e c) 2 ( 7 2 - l) = (m e c) 2 • 3 — ~4yv ( 28 ) 

K 2 5 j 



where (.) denotes the average over the distribution (27). Let W = -^^jy be the average particle 
energy. Then, the cooling rate due to synchrotron losses is 



fdW\ B 2 K 3[h) -<r 

H?- ^-ca T —e >^4-e"M (29) 

where the exponential is an approximate expression to include the suppression of cooling due to 
the Razin-Tsytovich effect. (7) = K^(l/Q) / K2(l/Q) — is the mean Lorentz factor of the pairs 
(see eq. [38]) and 



1 
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™ K211 ° l^J [a) ^ (30> 

is the Razin-Tsytovich Lorentz factor (Crusius & Schlickeiscr 1988). The cooling rate due to 
inverse-Compton scattering of accretion disk photons in the Thomson limit is 



Vdr) lc = - meC is^ r {—) ^7iy J dRR@ °w x4 (3i) 
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where R is the distance of a point in the accretion disk to the centre, Ri n and R out are the 
boundary radii of the disk, ®d(R) is the dimensionless temperature of the accretion disk material 
at radius R, and x = \/R 2 + z 2 . If the point source approximation (z 3> Rout) holds, Eq. (31) 
reduces to 



VdFsu- 



^{l + (5 v fz{tf K2 (iy 



(32) 



Fig. 7 illustrates that for a typical Shakura-Sunyaev disk of total luminosity L = 10 46 erg s 1 , 
surrounding a black hole of mass M = 1O 8 M , the point source approximation is an appropriate 
choice for z > 0.1 pc. The temperature (0 = 2 in Fig. 7) enters both Eqs. (31) and (32) 
in the same way and does therefore not have any influence on the validity of the point source 
approximation. 

The cooling rate due to pair bremsstrahlung emission of a thermal pair plasma has been 
calculated by Haug (1985b). We interpolate between his expressions for the extreme-relativistic 
and the non-relativistic regime: 



128 J„Jq 



for e < 1, 

96ar 2 m e c 3 nG(ln(2e) + 0.673) for 9 > 1. 
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For synchrotron-self-Compton scattering, we approximate the synchrotron radiation field to 
be distributed over a small range of energies and use its energy density 
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leading to 
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Using Eq. (8) for the adiabatic cooling rate, we find 
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The heating rate due to stochastic interactions with Alfven waves, according to Eq. (12), is 
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where v = (q + 3)/2. 

The mean kinetic energy of the electrons and positrons in a thermal plasma of temperature 

G is 
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(which in the non-relativistic limit approaches Wth = m e c 2 
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and the respective cooling rates are determined by 
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). Thus, 
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The annihilation rate in a thermal plasma is well approximated by Eq. (20). We found that 
numerical solutions of the Fokker-Planck equation in the case of a thermal pair plasma are in 
good agreement with the calculation described in this section as long as the pair temperature 
is > 4. At lower temperatures cooling of the thermal plasma is slower than predicted by our 
Fokker-Planck treatment. We attribute this mainly to the ultrarelativistic approximation used 
for the pair bremsstrahlung process which is not an appropriate approach for mildly relativistic 
temperatures. 

In Fig. 8, we show three examples of the evolution of a thermal pair plasma in an AGN 
jet with different initial densities. Plasmas of very high density end up with a density evolution 
which is almost independent of density, since initially the density evolution is dominated by pair 
annihilation losses which are proportional to the square of the density. At lower densities, the 
dilution of the plasma due to the expansion of the jet always dominates the decline of the density, 
which thus always remains well below the density of a plasma injected with high initial density. 
The temperature evolution depends only weakly on the density because it is governed by adiabatic 
losses. The differences in the temperature evolution curves are a consequence of the density 
dependence of bremsstrahlung losses and of the heating rate via v 2 a oc rC 1 . 

The x-ray and 7-ray emission of a dense, mildly relativistic pair plasma in the outer parts of 
an AGN jet is calculated in the same way as for the dilute thermalized pair plasma, as described 
at the end of the previous section. 
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5. Application to PKS 0208-512 

The quasar PKS 0208-512 is the most prominent example of MeV blazars, i. e. blazars 
exhibiting variable emission around several MeV and reaching extremely high fluxes in this 
frequency range, often dominating the total energy output of the source. Analyzing COMPTEL 
data of Phase I and II, Blom et al. (1995) found that during Phase I, PKS 0208-512 did not show 
a significant signal in the COMPTEL energy range, while in Phase II, a pronounced bump at 1 
- 3 MeV was detected. Quasi-simultaneous EGRET observations revealed a strong signal from 
PKS 0208-512 during both Phase I and II. However, in Phase II, EGRET detected the source 
with lower significance than during Phase I. The EGRET spectra appear to be harder during 
Phase I than in Phase II. During Phase III PKS 0208-512 is again detected by EGRET at a high 
flux level, whereas the COMPTEL data for that period do not show evidence for emission from 
PKS 0208-512 (Blom 1996). The variability in the EGRET range correlates with spectral index, 
which itself varies more strongly than in the case of the average blazar (see Stacy et al. 1996 
and Pohl et al. 1997). In fact, the flux level at 100 MeV stays approximately constant while the 
emission above 300 MeV is highly variable. 

The variability timescale of the high-energy component of the 7-ray spectrum of our model 
system is determined by the light travel time through the plasma component and corresponds to 
~ 6 h in the observer's frame. The timescales of quasi-thermalization and of dilution of the pair 
plasma due to pair annihilation and the expansion of the jet are of the order of a few weeks. An 
outburst at MeV energies due to transrelativistic, quasi-thermalizing pair plasma in a relativistic 
jet is therefore generally not expected to be correlated with a 7-ray outburst observed by EGRET 
within the same observing period. In general terms, the flux in the EGRET range reflects the 
injection rate of high-energy pairs, whereas the flux of X-ray and soft 7-ray photons is a measure 
of the injection rate of pairs at all energies averaged over several weeks. 

The instantaneous pair annihilation and pair bremsstrahlung radiation of one single 
thermalizing pair plasma blob is too weak to produce the observed MeV bump in the spectrum of 
PKS 0208-512. Therefore, if the MeV blazar phenomenon is to be explained by pair annihilation 
radiation of cooled, thermalized pair plasma in a relativistic jet, this jet must be quasi-continuously 
filled with pair plasma. 

If, in contrast, the MeV bump is produced by inverse Compton scattering of external radiation 
by the quasi-thermal, transrelativistic pair plasma in the jet, reacceleration of particles must be 
efficient because else the X-ray and soft 7-ray spectrum would just be the extension of the 7-ray 
power-law, resulting from Comptonization of soft radiation in the Compton-cooled nonthermal 
plasma. The bulk Compton process is therefore mainly relevant in the low to medium density case, 
in which relativistic temperatures can be maintained. However, a problem with this interpretation 
is that the quasi-equilibrium temperature of the pair plasma depends extremely sensitively on 
the magnetic field strength, the level of hydromagnetic turbulence and the injection height of 
the plasma blob, as indicated by Eq. (19). Therefore, if inverse-Compton scattering of external 
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radiation by quasi-thermal pair plasma in the jet plays an important role, we would expect to see 
bulk-Compton bumps in different objects at a largely different frequencies, for which no evidence 
has yet been found. There is no obvious reason why this bump should be produced preferrably at 
several MeV. 

The dashed line in Fig. 9 shows a fit to the EGRET spectrum during Phase II (VP 220) 
where we adopted the following set of parameters: n e = 5 • 10 5 cm -3 , Rb = 6 • 10 15 cm, ji± = 500, 
72± = 2 • 10 4 , s = 2.4, B = 0.7 G, b = 2, z, t = 2 • 10~ 3 pc, T = 10, 9 obs = 6°, L = lO^ergs -1 , 
Mq = 10 8 Mq. The high-energy spectrum results predominantly from inverse-Compton scattering 
of external (accretion-disk) photons. The photon spectrum has been integrated over the cooling 
time of the ultrarelativistic pair population because this time is only several hundred seconds 
(in the observer's frame) which cannot be resolved by COMPTEL and EGRET. The fluence 
resulting from our time integration has been re-converted into an average flux by multipling by 
the repetition rate of subsequent blobs. Since during an observed time interval At b s the blob 
travels a distance Al = DT [3t c Ai b s cm, a jet filled quasi-continuously with relativistic material 
(Al ps Rb) corresponds to a repetition time in the observer's frame of Ai b s ps 2000 s which 
translates to a quasi-continuous energy input into the jet of Lj et ps 5 • 10 43 erg s _1 during an 
EGRET outburst. 

The bump of the radiation from an accretion disk of L = 10 46 erg s _1 at redshift z ps 1 
corresponds to v F v ps 10~ 6 MeV cm~ 2 s _1 , if we assume Ho = 75km/(s ■ Mpc), qo = 0.5. No 
simultaneous broadband measurements in the radio to UV frequency range are available. For 
this reason, we make no attempt to fit the synchrotron component. Anyway, the radio emission 
probably originates in a much larger volume than considered here, in accord with the fact that 
there is no correlation between the simultaneously observed gamma-ray flux and cm-radio flux of 
blazars (Miicke et al. 1996, 1997). 

Using the methods described in the previous sections, we follow the evolution of the 
plasma blob through the transrelativistic phase until it is diluted by pair annihilation and the 
expansion of the jet to less than 10~ 4 of its initial density. Due to the rather low initial density, 
quasi-thermalization occurs at a temperature of ps 100, implying that the X-ray and soft 7-ray 
spectrum is strongly dominated by inverse Compton scattering of accretion disk radiation. 

In order to reproduce the MeV bump, we assume that some fraction of the jet is filled with 
pair plasma, evolving in the same way as the ultrarelativistic component, over a long period 
of time. We find this fraction to be ~ 10 %. This is consistent with the observation of strong 
variability in the EGRET regime if one assumes that during the quiescent phases the injection of 
relativistic pair plasma into the jet is very inefficient, and is in perfect agreement with the typical 
flaring duty cycle of EGRET detected blazars. The combined inverse-Compton, pair annihilation 
and pair bremsstrahlung spectrum emitted by the quasi-thermalized pair plasma jet as it moves 
out and is diluted due to jet expansion and pair annihilation, is shown by the solid line in Fig. 9. 

We point out that the MeV bump can be fit with a self-consistent pair annihilation spectrum 
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only under the assumption of a cylindrical jet, which seems to contradict the assumption that the 
gas pressure of the pair plasma greatly exceeds the magnetic pressure, most probably responsible 
for the jet confinement. 

6. Summary and conclusions 

Motivated by the standard model for 7-ray emission from blazars invoking jets filled with 
ultrarelativistic pair plasma oriented at a small angle with respect to the line of sight, we 
investigated the further evolution of the leptonic material inside such a jet after the initial phase 
in which the high-energy (> 100 MeV) 7-ray radiation is produced. 

We demonstrated that either the action of elastic scattering of particles off each other or the 
balance of reacceleration by turbulent plasma waves to radiative losses can, under a wide range of 
parameters, establish a quasi-thermal particle distribution inside the jet, irrespective of the initial 
pair distribution at the time of injection. Extreme conditions (high density, very weak magnetic 
field, injection close to the accretion disk) are necessary to achieve mildly relativistic temperatures 
(6 ~ 1). Radiation from a quasi-thermal pair plasma can produce the temporary MeV bumps 
observed in MeV blazars. This can be accomplished most easily via inverse-Compton scattering of 
external soft radiation by a quasi-continuously filled jet. 

Alternatively, the MeV bump could be produced by pair annihilation radiation, if cooling 
of the pair distribution is very efficient and the dilution of the plasma through jet expansion is 
negligible. We found that only the unrealistic assumption of a cylindrical jet yields an acceptable 
fit to the observed MeV bump, if the pair plasma producing this bump is subject to the same 
processes (in the same environment) as the ultrarelativistic pair plasma producing the EGRET 
spectrum. 

In this picture, an MeV outburst reflects an increasing supply of relativistic electrons and 
positrons into the jet on much longer timescales than those typically observed in EGRET outbursts 
at higher energies. It is only weakly or not at all correlated to activity in the EGRET range which 
corresponds to an increased energy density at the acceleration site, leading to a harder particle 
spectrum of the injected pairs. 

If the bulk Compton process is the dominant mechanism for the production of the MeV bump 
in MeV blazars, then there is no obvious reason for a "universality" of the bump photon energy 
at 3 — 10 MeV. Mainly depending on the magnetic field strength, the level of hydromagnetic 
turbulence and the injection height of the plasma blob above the accretion disk, similar sources 
with bulk-Compton radiation bumps at different energies should exist as well. 

In general, our treatment is also applicable to the so-called hadronic jet models, where 
ultrarelativistic protons initiate a pair cascade in the jet (Mannheim et al. 1991, Mannheim & 
Biermann 1992, Mannheim 1993). However, although the jet simulation code of BMS is able to 
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handle problems in which the blob is optically thick to 77 pair production, it is clearly beyond the 
scope of this paper to follow the proton-initiated pair cascades in detail. Nevertheless, a few general 
conclusions may be drawn from the fact that in those models, the magnetic fields are fairly high 
and a considerable level of turbulence might be present. Our simulations indicate that at magnetic 
fields of B ^ 10 G, as usually assumed in hadronic jet models, and the moderate pair densities 
expected to result in the cascades, the pair plasma will maintain a highly relativistic temperature, 
^ 100. Thus, no pair annihilation feature will be observable; instead, strong synchrotron and 
SSC bumps at ~ 3 ■ 10 7 -D/(l + z) (B/G) O 2 Hz and ~ 4 • 10 8 -D/(l + z) (B/G) 4 Hz, respectively, 
would be expected from the thermalized cascade plasma. 

We thank the anonymous referee for valuable comments and suggestions which led to 
significant improvements of the manuscript. This work has been partially supported by NASA 
grant NAG 5-4055. R. S. acknowledges partial support by the DARA. 

A. Energy exchange and dispersion rates 
A.l. Energy dispersion due to elastic scattering 

Following Dermer (1985), the energy dispersion rate is given by 
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Here 7,. is the Lorentz factor of the relative motion of the interacting particles, 7 C is the Lorentz 
factor (normalized velocity: (3 C ) of motion of the center-of-momentum (cm) frame with respect 
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to the blob frame, 7 cm is the Lorentz factor (normalized velocity: (5 cm ) of the particles in the 
center-of- momentum frame, and u is the cm angle cosine between f3 cm and (5 C . 

It is accurate within a few percent to retain only the leading terms ~ u 2 and ~ (1 — u 2 ), 
respectively, that is the terms ~ u 2 ^ 2 ./2 and ~ (1 - u 2 )j 2 (ln\/2/2 + lnA)/2 in Eq. (6). 
Nayakshin and Melia (1998) suggested an approximation based on the further neglect of all 
terms not containing an In A, arguing that this contribution is the one resulting from small-angle 
scatterings and that, at the point where other contributions may become important, the use of the 
Fokker-Planck equation is no longer a good approach since then large-angle scattering, changing 
the particles' energy considerably, dominate the energy exchange in elastic scattering events. Our 
simulations show that in the case of the system we are dealing with the energy distributions of 
the pairs remain relatively narrow. In this case, the dispersion may indeed be approximated 
by the expression found by Nayakshin & Melia which ensures the validity of the Fokker-Planck 
equation. Nevertheless, we retain both contributions (~ u 2 and ~ [1 — u 2 }) in order to see whether 
large-angle scattering is important at the point where the plasma starts to thermalize (i. e. when 
elastic scattering becomes the dominant process). 

With this simplification, the dispersion rate due to elastic scattering off an isotropic 
distribution of particles 712(72) is given by 
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is the monoenergetic dispersion rate and 
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A. 2. Bremsstrahlung energy loss 

Using again the technique of Dernier (1985), we find the single-particle energy loss rate due 
to pair bremsstrahlung emission to be 
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(cf. DL). The maximum photon energy e™ x is given by 
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A. 3. Energy dispersion due to inverse-Compton scattering 

The energy dispersion rate due to inverse-Compton scattering of photons of differential 
density n p h(e, £l p h) in the Thomson regime is 
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where rj p h is the cosine of the angle between soft photon momentum and jet axis. Inserting the 
photon number density (1.8) yields 
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where Ri/ a is the radius of the inner/outer edge of the accregion disk, respectively, x = V R 2 + z 2 
andC(n)=Er=ifc-". 

In the case of an isotropic soft photon field we have n p h(e, Q p h) = n p j l (e)/Ai: (e. g. synchrotron 
photons in our model), and Eq. (19) reduces to 
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A. 4. Pair annihilation losses 

The general expression for the catastrophic particle losses due to pair annihilation can be 
written as 
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is the angle-averaged reaction rate and 
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(Svensson 1982). 7™™, l^ x are given by the kinematic restrictions, namely 
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B. Approximative analytical solution 

As we found in section 4, for ultrarelativistic particles with 7 > 100, and neglecting pair 
annihilation, the Focker-Planck equation (1) can be approximated as 
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where n = 71(7, i) and <9n/(x[i]<9t) = dn/dT. Laplace transformation of Eq. (Bl) to 
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In the first step we search for a homogeneous solution ^.(7) which satisfies the equation 

{phW)' -g( 7 )h = 0. (Bll) 

Substituting 

h{ 1 )=p{ 1 r l ' 2 P{x), (B12) 
choosing x = 7 _1 and inserting the functions g and p, yields a differential equation for P(x): 

x P"(x) + 2P'(x)-{(x + rj) P(x) = (B13) 
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Two linearly independent solutions of Eq. (B 13) are given by 
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where a = y/Q. Since we did not find an analytical solution for the inverse Laplace transformation 
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of the Green's function resulting from the homogeneous solutions ^1/2(7) = P^ l ^ 2 {l) -^1/2(1/7) 
we simplify Eq. (B 13) in a convenient manner. Since the energy loss and dispersion rates are 
dominated by the inverse-Compton process, we find for the coefficients 
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where e and n p h{e) are dimensionless energy and spectral photon number of the accretion disk 
photon field in the blob rest frame. Hence, we can estimate 

~ -L ~ 10 5 (B21) 
A 2 (e) 

for the accretion disk model we used (see BMS, section 2.1). Since for the Laplace transform 

> 0, we have |£| >, 10 10 and rj ~ 10 5 . Thus, for particle energies 7 < 10 4 , (x 3> 77, and we may 
neglect the term ~ r\. The numerical values for the coefficients A\ and A 2 which we find during 
our simulations confirm the above estimate. Therefore, Eq. (B 13) is well approximated by 

P"(x) + ^P'(x) - C P{x) = (B22) 

which is solved by 

P l/2 {x) = \e ±ax (B23) 
where, again a = The homogeneous solutions are thus 

h 1/2 h) = e(^ ±a )^- 3 . (B24) 
The Laplace-transformed Green's function G(j,jo,s) is then constructed as 
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where 
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and W{hi,h 2 } denotes the Wronskian. The functions h\ and h 2 are chosen in a way that h\ 
fullfills the boundary condition at 7 — > 1 and h 2 fullfills the boundary condition for 7 — > 00 where 
we simply have to use the condition of finite solutions G(7, 70, s)//(7o) at 7 ^ cx) and 7^0. 
(In the approach used here we can also write Eq. (Bl) differential in kinetic energy E := 7 — 1 
yielding the same solution in E instead of 7; then the boundary condition is determined by a finite 
value of G/f at E = 0.) Thus, the resulting Laplace transformed Green's function is 
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The inverse Laplace transform of this Green's function is 
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and an approximative solution of Eq. (B 1) is given by 
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Fig. 1. — Energy-loss rates in a dense relativistic jet for the following parameters: n = 3 • 10 6 cm 
(7) « 30, B = 0.1 G, 5B/B = 0.1, z = 10" 3 pc, L = lO^ergs" 1 , T = 15 
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Fig. 2. — Energy-dispersion rates in a dense relativistic jet for the same set of parameters as in 
Fig. 1 
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Fig. 3. — Evolution of the pair distribution functions through the transrelativistic phase in the 
case of high particle density, uq = 10 9 cm~ 3 . Initial conditions: ji± = 200, 72± = 2 • 10 4 , s = 2, 
R B = 5- 10 12 cm, B = IG, Zi = 10~ 3 pc; L = 10 46 erg s _1 , M = 10 s M Q ,T = 10. The solid line 
indicates the state of the plasma after the ultrarelativistic treatment described in BMS 
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Fig. 4. — Comparison of the final state of the simulation shown in Fig. 3 (solid) to a thermal 
spectrum of temperature = 0.83 (dotted) 
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Fig. 5. — Evolution of the pair distribution functions through the transrelativistic phase in the case 
of moderate particle density, no = 10 5 cm -3 . Initial conditions: ji± = 500, j2± = 3 • 10 4 , s = 2.1, 
R B = 2- 10 15 cm, B = IG, Zi = 10~ 3 pc; L = 10 46 erg s -1 , M = 10 8 M , T = 10. The solid line 
indicates the state of the plasma after the ultrarelativistic treatment described in BMS 
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Fig. 6. — Comparison of the final state of the simulation shown in Fig. 5 (solid) to a thermal 
spectrum of temperature 6 = 1.3 (dotted) 



-35- 




z[pc] 



Fig. 7. — Cooling rates of a thermal pair plasma due to inverse-Compton scattering of accretion 
disk photons as a function of distance from the disk. Solid: Extended source; dashed: point source 
approximation. L = 10 46 erg s" 1 , M = 10 8 M , = 2 
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Fig. 8. — Evolution of thermal pair plasmas in a relativistic jet. Initial temperature 6o = 1; 
B = 0.1 G, z = 10" 2 pc, L = 10 46 erg s _1 , V = 10, 5B/B = 0.1, B(z) oc z -1 ; initial densities 
no = 10 s cm" 3 (dotted), 10 9 cm" 3 (solid), and 10 10 cm" 3 (dot-dashed), respectively. Thick curves 
show the temperature evolution, thin curves the density (multiplied by 10" 9 ) 
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Fig. 9.— The COMPTEL and EGRET spectrum of PKS 0208-512 during Phase II. Dashed: high- 
energy (predominantly external inverse-Compton) spectrum from ultrarelativistic plasma blobs; 
solid: inverse-Compton spectrum from the later stages of the jet evolution, under the assumption 
of ~ 10 % of the jet being filled with transrelativistic pair plasma 



